arXiv:1504.05989vl [cond-mat.stat-mech] 22 Apr 2015 


CCNY-HEP-15-01 


Statistics of two-dimensional random walks, 
the “cyclic sieving phenomenon” 
and the Hofstadter model 

Stefan Mashkevich0 

Schrodinger, 120 West 45th St., New York, NY 10036, USA 

and 

Bogolyubov Institute for Theoretical Physics, 03143 Kiev, Ukraine 

Stephane Ouvryi 

Physics Department, City College of the CUNY, New York, NY 10031, USA 

and 

Laboratoire de Physique Theorique et Modeles Statistiques, CNRS-Universite Paris Sud, 
Faculte des Sciences d’Orsay, 91405 Orsay, France 

Alexios Polychronakod^ 

Physics Department, City College of the CUNY, New York, NY 10031, USA 

and 

The Graduate Center, CUNY, New York, NY 10016 
April 24, 2015 

Abstract 

We focus on the algebraic area probability distribution of planar random walks on a 
square lattice with mi, m 2 , h and I 2 steps right, left, up and down. We aim, in particular, 

2i7r 

at the algebraic area generating function evaluated at Q = e « , a root of 

unity, when both mi—m 2 and I 1 —I 2 are multiples of q. In the simple case of staircase walks, 

2i7r 

a geometrical interpretation of ^m,o,«,o(e ^ ) in terms of the cyclic sieving phenomenon is 
illustrated. Then, an expression for which is relevant to the Stembridge’s 

case, is proposed. Finally, the related problem of evaluating the nth moments of the 
Hofstadter Hamiltonian in the commensurate case is addressed. 
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1 Introduction 


The cyclic sieving phenomenon [T] is quite ubiquitous in combinatorics, whereby some 
hnite sets have both a cyclic symmetry and a generating function which, when evaluated 
at roots of unity, happens to count some symmetry classes of the sets. A paramount 
example is the collection of the /-subsets of a (m -|- /)-set with the Q-binomial generating 
function 

+ l\ _ [m + /]q! 

I 


m 


Q-Wq- 


( 1 ) 


The Q-factorial is defined as 


[^]q' — n 1 _ n “ 1(1 -h Q)(1 Q Q^) • ■ ■ (1 Q 

1 ^ ^ 

2=1 


+ q'-') . 


( 2 ) 


Denoting by c the cycling generator (one-step cyclic permutation of the (m -|- 1) set), one 
hnds that evaluated at Q = counts the number of /-subsets fixed by (i.e., 

invariant with respect to) for p = 1, 2,..., m -|- /. In the particular case Q = —1 one 
refers to the Stembridge’s phenomenon |2]. The fact that inte^rs show up here stems 
from a well-known identity: for any integers p, q mutually prim^ and m a multiple of g, 
one has 

m-\-l 
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A planar random walk on a square lattice is dehned as an ordered sequence of steps 
to the right x, left up y, and down y~^, with numbers of steps mi, m 2 , /i and 
I 2 , respectively. We will refer to such a walk as an (mi, m 2 ,/i,/ 2 ) walk. The walk is 
closed if mi = m 2 = m and /i = I 2 = n/2 — m, where n is the total number of steps 
(which is necessarily even), and open otherwise. We extend the standard dehnition of 
the algebraic area enclosed by a closed walk onto an open walk, as follows: close an open 
walk by connecting its end point with its start point, adding on to the end of the walk 
the minimum necessary number of steps, first vertical, then horizontal. E.g., if mi > m 2 
and /i > / 2 , then we close the walk by adding /i — I 2 steps down followed by mi — m 2 
steps to the left. The area of the open walk is defined as the area of the closed walk thus 
obtained; essentially, it is the algebraic area under the open walk. 

The generating function of algebraic area probability distribution of 

(mi, m 2 ,/i,/ 2 ) walks originating from a given point on the lattice is defined in terms of 
the number of such walks with algebraic area A: 




( 4 ) 


A=—c 


■^In the sequel, Q denotes the argument of the Q-binomial, whereas the integer q denotes the root of 
unity Q = exp(2i7r/(7) at which the generating functions are evaluated. 
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Zmi,m 2 ,h,i 2 i.^) can be viewed as a generalization of the Q-binomial coefficient: If x and y 
(identified with steps to the right and up, respectively) satisfy the commutation relation 
xy = Qyx, then 

{x + y + x~^ + y~^)^ = Zm,,m 2 ,iui 2 {Q)y~^^y^^x'^^x~"^^ . (5) 

mi,m2,h,h 

mi+m2+h+l2=n 

It can be seen easily that ^mi,m 2 ,Zi,Z 2 (Q) satishes the recurrence relation 


,m2 


,ll,l2 (Q) —1,/2 (Q) ■^mi,m2,/i,/2 —1 (q) 




h—h ' 


( 6 ) 


with the initial condition ^o,o,o,o(q) = 1- If is therefore a polynomial in Q and Q“^, the 
solution of Eq. ([6]), but whose explicit expression remains unknown. 

Trivially, when evaluated at Q = 1, 


■^mi,m2,Zi,Z2 (I) 


mi + m2 + h + h 

mi,m2,h,l2 


(7) 


is the usual multinomial coefficient that counts the number of (mi, m 2 ,/i,/ 2 ) walks. 

Less trivially, for “biased walks”, i.e., ones that can go right, up and down, but not 
left — that is to say with m 2 = 0 — the generating function Zmi,o,q,i 2 (Q) been found 
[3] to be 


min(Zi,Z 2 ) r 


■^mi,0,Zi,Z2 (Q) — ^ ^ 


fc =0 


mi + /i + ^2^ /mi + /i + /2 
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mi + li — k\ f mi + I 2 — k 


mi 


mi 


( 8 ) 


Again, when evaluated at Q = 1, the multinomial counting for the number of (mi, 0, /i, I 2 ) 
biased walks is recovered: 


■^mi,0,Zi,Z2(l) — 


mi + /i + /2 


mi, /i, I 2 

When also h = 0, ^mi,o,Zi,o(Q) yields, as it should, the Q-binomial generating function 


(9) 


=(“ + ') (10) 

\ ^ Jq 

for the probability distribution of the algebraic area under staircase walks, ones that can 
go only m steps right and I steps up. One notes that this Q-binomial has been already 
introduced in Eq. ([1]) and evaluated at Q a root of unity in Eq. ([3]) in the context of the 
cyclic sieving phenomenon for the /-subsets of the (m -|- /)-set. 

These considerations evoke a natural question: Can Zmi,o,Zi,Z 2 (Q) Eq- dH]) with Q a 
root of unity be an integer, and if so, what does this integer count? More generally, for 


3 




{mi, m 2 ,1 1 , h) walks, despite ^mi,m 2 ,/i,« 2 (Q) being generally unknown, can it be evaluated 
at Q a root of unity; if so, can it yield, at least in certain cases, an integer; if yes, what 
does this integer count? 

Apart from these cyclic sieving combinatorics/counting considerations, ^mi,m 2 ,Zi,i 2 (Q) 
happens to be of interest not only for random walks but also for the quantum Hofstadter 
problem [1], thanks to the formal mapping [5] between the algebraic area generating 
function for closed random walks of length n 

nl2 

■^n(Q) ^ ^ ^m,m,^—m,^—m{Q') (H) 

m,=0 

and the n-moments of the Hofstadter Hamiltonian 

Z^{d^) = TiH:; , ( 12 ) 

where 7 = 27r0/0o is the flux per plaquette in units of the flux quantum and Q has been 
taken to be e*'*'. Of particular interest is the commensurate case 7 = 27rp/q, with p and 
q relative primes. Therefore, another motivation in evaluating when Q is 

a root of unity stems from the Hofstadter quantum spectrum itself for a commensurate 
flux. 

In this regard, for “closed” biased walks of length n — “closed” here being defined by 
h = h = {n — mi)/2 and the choice of particular boundary conditions on the horizontal 
axis — an explicit quantum mapping has indeed been found [ 6 ]: 


L q J 

= Tr 7 = 27r/g , 

m=0 


(13) 


where now stands for a truncated Hofstadter Hamiltonian with the horizontal hopping 
to the left absent. In the sum above, each .^gm,o,(n-gm)/ 2 ,(n-gm )/2 has to be understood 
with n and qm being of the same parity (otherwise it has to be taken equal to 0). Note 
also that the number of right steps qm is a multiple of q, so that one can view the “closed” 
biased walks as winding on a cylinder with circumference q. 


2 The sieving phenomenon and staircase random walks 

As a warm-up, let us rephrase the staircase random walk algebraic area generating function 
^m,o,i,o(Q) evaluated at Q a root of unity in the context of the cyclic sieving phenomenon. 

One first notes that an (m, 0, 1, 0) staircase walk consisting of m steps to the right and 
I steps up uniquely corresponds to an Z-subset of the set {1,... ,m -f /}. For example, 
xxyxxy, which is one of the possible (4, 0, 2, 0) walks, corresponds to the subset {3, 6} of 
the set {1,..., 6 }. A cyclic permutation c of the above-mentioned set is 1 —)■ 2 —)■ 3 —)■ 

2i7rp 

...—)■ m -(- Z ^ 1. The cyclic sieving phenomenon is the statement that Zm,o,i,o{e'^+‘) 
equals the number of the subsets fixed by for p = 1, 2,..., m -f Z. 
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Next, one notes that when both m and I are mnltiples of q, Eq. (|3]) implies 


■‘m,0,1,0 


(e 




( 14 ) 


and the qnantity on the RHS connts the number of staircase walks with m/q steps to the 
right and Ijq steps up. 

This counting allows for a simple interpretation in terms of the staircase walk (m, 0, /, 0). 
The most general /-element subset of the set {1,... ,m -|- /} hxed by d’ has the form 

{/l, *2, ■ ■ ■ , b/q, +P, *2 +P, • • • , b/g +P, ■ ■ ■ , *1 + (? - l)p, ^2 + (? “ 1)P, ■ ■ ■ , H/q + (? “ l)p}, 

where q= {m + l)/p (i.e., q copies of a subset {R, *2, • • •, b/g}, with b/g < P, shifted by p 
with respect to each other). A subset {/i,/2, • • • di/q} corresponds, as formulated above, 
to a (^,0,^,0) walk. Thus, any (m,0,/,0) walk corresponding to a subset hxed by d 
is a repetition of q identical “building blocks”, shifted with respect to each other by p 
steps. Equation (ITT)) expresses the fact that the total number of such walks is equal to 
the number of possible building blocks. 

Consider an example: m = 9, / = 3. The values of .^ 9 ,o, 3 ,o(e^) with p ranging from 1 
to 12 are: {0, 0, 0, 4, 0, 0, 0, 4, 0, 0, 0, 220}. The hrst 4 in this sequence, corresponding 
to p = 4, counts the number of ways that a (9, 0, 3, 0) walk can be constructed from 3 
building blocks, each block being a (3, 0,1, 0) walk: there are 4 such blocks possible. Each 
of those walks is hxed by c^. The second 4, corresponding to p = 8 , rehects the fact that 
any subset hxed by is also, trivially, hxed by c®. Finally, every walk is hxed by = I, 
so the 220 is the total number of walks ( 3 ). The zero values of ) rehect the 

fact that there are no walks hxed by c, c^, c^, c®, etc. 

2i7rp 

Note also that the non-zero values can be directly retrieved from Zm,o,i,o{G~) with 
p = 1,2, ... ,q where q is the GCD of m and /. In the example above .^ 9 ,o, 3 ,o(e“^) with 
p = 1, 2, 3 directly yields {4, 4, 220}. 


3 The Stembridge phenomenon for biased and unbi¬ 
ased random walks 

Consider now in Eq. ([H]) the generating function Zmi,o,h,i 2 {Q) for fhe algebraic area dis¬ 
tribution of biased walks and evaluate it with Q a root of unity. Looking either at the 
staircase walks counting in Eq. ([3]) with m a multiple of q or at the g-periodic sum in 

2i7r 

Eq. ffl3]) . let mi be a multiple of q. It follows immediately that ^mi,o,q,Z 2 (o ■? ) is an inte¬ 
ger. Furthermore, as for staircase walk^, all Zmi,o,Zi,Z 2 (o~)’s for p = 1, 2,..., g are also 
non-vanishing integers. 

One can go a step further by having not only mi but also |/i — / 2 I be a multiple of g, 

2i7rp 

®The analogy stops here: Zm^,o,h,h (e™i+'i+*2) with p = 1, 2 ,..., mi +li +I2 would yield non-integer 
values for certain values of p. 
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as suggested by Eq. (fH]) . Some algebra allows then to rewrite 



as 


k) 

( m\+h—k \ / mi+l2—k \ 

’f )\ ^ )' 

In Eq. ([15]), Efc=min(Zi means that the sum over k is from min(/i,Z 2 ) down to 0 by 
steps of minus q — from which one deduces that the entries of the last two binomials are 
indeed integers. It is also understood that when mi = 0, the ratio is 

non-vanishing and actually equal to 1 only when k = min(/i, I 2 ). This yields, as it should, 

~ number of {li^h) walks on the vertical axis which 

2i7r 

all have, trivially, a vanishing algebraic area so that « ) = Zo,o,Zi,Z 2 (l)- 

In 0151) for eachl§ k the first binomial counts the number of ways to pick k elements 
out of the (mi -|- /i -I- / 2 )-element set, whereas the last two binomials are the numbers of 
(^, and (^, staircase walks, respectively. Finally, the ratio in fflSj) encodes 
the fact that the li — k steps up and I 2 — k steps down “annilihate” each other on the 
vertical axis so that these staircase walks cannot be considered as independent. 

Now, let us focus on the Stembridge’s case Q = —1, that is to say q = 2: one can then 
simplify (1T5|) to get the quite convincing expression 


^mi,0,Zi,Z2 


^ 2i7r 

e 


) = 


0,-9 

E 


'^1 + + ^2 
k 


mi (mi -|- Zi -|- Z 2 — 2,k) 
{mi + li- k){mi + I 2 - 


^mi,0,Zi,Z2 


(- 1 ) 




mi+h+l2 

2 

Z 1 +Z 2 

2 


(16) 


Furthermore, one can go a step further and relax the biased constraint m 2 = 0, that is 
to say consider now more general (mi, m 2 , h, h) walks, with both |mi — m 2 | and |Zi — / 2 I 
multiples of q. Again, in the Stembridge’s case q = 2, one obtains 


■‘mi,m2,h,l2 ) 


- 1 ) = 


mi -|- m 2 
mi 


li + I 2 
h 


mi+m2+Zi+Z2 

2 

Z 1 +Z 2 

2 


( 17 ) 


which does reduce to Eq. (ITB]) for m 2 = 0 biased walks. 

Eq. (fT7|) . and consequently Eq. (IT^ . yields integers, since |mi — m 2 | and |/i — / 2 I 
being even implies that mi -|-m 2 and h + I 2 are even as well. A Stembridge combinatorial 

®Note that for each such k, 

mi{mi+h+h-2k) 

{mi + h - k){mi + h - k)\ ir^ J \ J 

is by itself an integer, which is nothing but saying that, trivially, 

a(a -1-6-1-c) /a-|-6\/a-|-c\ 

(a -I- b){a -l-c)V a / y a ) 

is an integer for any integers a, 6, c. 
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□□ 
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□□ 
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□□ 

010000 

□□ 


^1 + ^ 2 + ^1 + ^2 
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right 


0 = up 

0 = down 


Figure 1: A Stembridge interpretation. 


interpretation of these integers is depicted in Fig. 1. Divide the {mi+m 2 + li + l 2 ) cells into 
two equal subsets (g = 2, i.e., two building blocks). Then consider four types of objects 
(corresponding to four directions — up, down, left, right) and dehne a subset as “hxed in 
the weak sense” by c(™'i+™' 2 +^i+i 2)/2 — meaning that when acting by c, “up” is considered 
to be the same as “down” and “left” the same as “right”. The formula is then interpreted 
as follows: the rightmost factor counts the number of ways that half of the (up + down) 
objects, (/i + ^ 2 ) 72 , can be distributed between half the cells, (mi + m 2 + + / 2 )/ 2 . Then 

two other factors distinguish between “up” and “down” (/i out of li + I 2 ) and between 
“left” and “right” (mi out of mi + m 2 ). 

The simple expression of Eq. flTT)) asks for a generalization to 

2i7r 

walks with both |mi — m 2 | and |Zi — / 2 I multiples of q. 

2 \-k 

Such an expression for Zmi,m 2 ,h,i 2 \^ ) would in turn allow one to retrieve the n-moments 

of the Hofstadter Hamiltonian at commensurate flux by summing ^m,m,/=n/ 2 -m,i=n/ 2 -m(Q) 
over m, like in Eq. flTT|) — with in that case both |mi — m 2 | and |/i — / 2 I vanishing, in 
principle an even more simple situation — and then using Eq. (IT^ . Following this line of 
reasoning for q = 2, i.e., the Q = —1 Stembridge case, one eventually gets from Eq. (Ell 
the nth moment of 


TtH: 


n/2 


n/2 


^ ^ ^m,m,l=n/2—m,l=n/2—m(^ 1 ) ^ ^ 


m=0 


m=0 



/2(t-m)Wt\ 
V I - m / \m J 


( 18 ) 


This last result is a sum rule for the Hofstadter spectrum at commensurate flux 7 = vr. 
Below we compute TrFf^Tr/g particular values of g = 2, 3,4, 6 directly from the quantum 
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mechanics formulation (i.e., without relying on an evaluation of '^ ))• 

4 The nth moments of the Hofstadter Hamiltonian 

Up to now we have used the mapping of the quantum Hofstadter problem on classical 
random walks to arrive at the Hofstadter sum rule (fT 8 |) . This sum rule can be directly 
retrieved from the quantum Hofstadter Hamiltonian itself, as can be other sum rules for 
simple values of q = 3,4,6. 

For concreteness, we wish to compute the nth moment 

Tr n” (19) 

where is the Hofstadter Hamiltonian 

= X + x~^ + y + y~^ , xy = e^'^yx ; (20) 

X and y stand now for the quantum lattice hopping operators on the horizontal and vertical 
directions respectively: they do non commute, due to the presence of the perpendicular 
magnetic held with hux 7 per unit cell. The dehnition of trace in flT^ is such that 

Tr(a;™|/") = 5^,o<5n.o. ( 21 ) 

We dehne the generating function 

00 

= = —= ( 22 ) 

n=0 ' ' 

where the last rewriting follows from the vanishing of Tr if” for odd n. 

From now on we focus on the commensurate case 7 = 27r/g, for which the operators 
x'^ and y^ do commute with x and y and are Casimirs, and so is 

Cq = x‘^ -\- x~^ + y^ + y~^ ■ (23) 

Cq is, essentially, the Hofstadter Hamiltonian with q = I, i.e., 7 = 27r, that is to say zero 
hux. One easily obtains 

TrCf =((")' (24) 

and 


Tr HL/, = 4 , 

Tr iiaTT/q = 28 + 8 cos(27r/g) , 

TrUf+' = Trii2;+/ = 0, (25) 

so Tr will be nonzero only for even values of n, as expected. 

The operators x, y can be represented as q x q matrices (with x'^ and y'^ proportional 
to the identity matrix) and therefore so can H 2 n/q- As a result, it satishes a characteristic 




equation of degree q, and so can be expressed in terms of lower powers of H 2 n/q 

and the Casimir Cq. Defining the “parity” of the monomial as (—I)™'"*'”, has 

parity (— 1 )"^ and only like-parity powers of -ff 27 r/<j will appear in the expression of 
Overall we have 

n/, = C, + 2,Hit + + ■ ■ ■ + " • (26) 

where are g-dependent numerical coefficients. Looking back at random walks, this is 
a rewriting of Eq. (jS]) but with the different terms regrouped into powers of H 2 n/q itself, 
which is possible for a rational flux 7 = 27r/g. The term Cq arises out of the terms 
x~'^, y~^ in the expansion of H^^^iq- The fact that the coefficient Cq -2 of the term 

is 2q can be proven using combinatorics and the identity 1 -|- -!-••• + = 0 . 

Performing the above rewriting explicitly for the hrst few values of q, we obtain 


g = 2: Hi =4 + 6*2 
q = ^ H ^^^3 = 6iL27r/3 + C's 

g = 4: Hi/, =shI/, + C,-4 

2n 

? = 5 : hf |^/5 = - (5 - 2 cos y) 6 ^ 2+5 + C'g 

g = 6 : Hl/^= 12Ht/, - 2AHI/, + A + C, (27) 

The above relations will allow for the exact evaluation of Tr The zero flux case 

g = 1, i.e., 7 = 27r (corresponding to the unweighted random walk, Q = 1) is explicitly 
solved as ^ 

TV . 2(1, t) = -K(U) (28) 

\n J TT 

with K{k) the complete elliptic integral of the hrst kind. 

The simplest nontrivial hux is g = 2, i.e., 7 = tt, that is, in the language of random 
walks, the Stembridge case Q = — 1. It serves to demonstrate a common phenomenon, the 
appearance of various different-looking expressions for the same quantity related through 
combinatorial identities. For the generating function we have 


= Tr 


1 - tmi 
1 


= Tr 


:Tr 


1 — f^(4 + 6 * 2 ) 
1 2 1 




:K 


TT 1 — 4f2 \1 — 




(29) 


where we have used the fact that Cq is essentially the free Hamiltonian H 2 -k- Equivalently, 


Tr Hl^ = 

TT 



(30) 
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An alternative formula is obtained by writing as a sum of two commuting parts, 
Hi = A +€2 = ^ ++ x~‘^ + 2 /^ + = (x + x~^y + {y + y~^Y (31) 

and thus 

Tri/f = Tr [{x + x-^f + {y + y-^]" 

= ^ W Tr (x + x-^)2^ Tr {y + 

k=0 ^ ' 

(n\ (2k\ (2n — 2k\ , , 


recovering Eq. dH. 

Yet another formula is obtained by noticing that the value of the Casimir C 2 is between 
—4 and 4 (since x^ and y"^ are phases) and thus it makes sense to put 


Co 


COSffl = 


(33) 


with some angle (/), it terms of which 


Hi = 4(1 + COS0) = 8cos^ — = 2 (e *2 + e *2 


(34) 


In this representation 


2n 


H2n ^ 2 - I + e-*2 ) = 


2"Y 


2n 


m 




2” ^ ^ j ~ > (35) 

m '' 


where in the last step we used the fact that only the real part of e®C "*)</> contributes. 
Changing the variable m = n + k, —n < k < n, using the formula 


cos k(j) = 2^ 




1=0 


k-l\ I 


(cos0)^-2' , fc>l 


(36) 


and using Eq. fl33|) . we obtain 


k=l 1=0 ^ 


I 


^k—2l 


(37) 


(in the above, is dehned to be 1 when both numerator and denominator are zero, as 
well as the ratios or below.) Taking the trace and using Eq. fl24)) . only even k 
values will contribute, and we obain a third formula for q = 2: 


[n/2\ k 

= EY(-1)'2 

k=^ 1=0 


n—2k-\-2l 


2n 


n + 2k J 2k — I \ I 


2k /2k - /2k - 21 


k-l 


(38) 
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It is remarkable that the above three different-looking formulas all agree and give the 
same result. The last formula may look a bit unmotivated and of dubious value, as it is a 
double sum (rather than a simple sum like the first two), but it is given in order to make 
contact with the upcoming results for higher q. 

For the case q = 3, the derivation proceeds by solving the characteristic equation for 
hf 27 r/ 3 , = 6 A -|- (Fa. The solutions are 


A„ = 2 V 2 


cos ■ 


0 -|- 2p7r 


p = 0,l,2 


(39) 


with the angle 0 now satisfying 


cos 0 


^3 

4 V 2 ' 


(40) 


Evaluating Tr amounts to averaging A” over the three values of p. In a calculation 

analogous to the third formula for q = 2 we have 


Tr 


p=0,l,2 P=0,l,2 


Tr 


<^-\-2p7T 


+ e~ 


■ 4>+2p7r 


2n 


3 ^ 

p,m 


2 nA 
m J 


0+2p7r 


rjr gi2(n-m)^ 


(41) 


The sum over p above will give zero unl ess n — m is a multiple of 3. Putting n — m = 3k, 
we have 


Tr = 2” 


2n 
n + 3k 


Tr = 2” ^ 


2n 

n + 3k 


Tr cos(2A:0) . 


(42) 


The rest of the calculation proceeds in a way similar to the derivation of the last formula 
for g = 2. The result is 


TiH 


2n 

27r/3 


[n/3] k 

ED 

A:=0 1=0 


- 1)'2 


I r)n—3k+3l 


2n 

n + 3k 


2k 

2k - I 



f2k-2lV 

V k-l ) 


(43) 


The similarity with fl38|) for g = 2 is obvious. 

For g = 4 we proceed in a similar way. The solutions to the characteristic equation 
A"^ = 8 A^ + 4 — C 4 are 


with 


Ap = ±Y 4 ± y /12 - 1 - C 4 


= 2 V 2 


cos 


0 -f 2p7r 


p = 0,1, 2, 3 


cos (j) 


4 + ^4 
8 


(44) 

(45) 
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Expressing Tr “ i Sp choosing either the hrst or the second expression 

above for Ap leads to two different-looking formulas. The hrst is 


[n/2] [k/2] 




2n—2k—4l ^k—2l 


k=0 1=0 


n 


k\ {21 


2k \2l \ I 


(46) 


The second one proceeds, again, in a way similar to the last of q = 2 and the q = 3 cases 
and gives 


[n/2] [k/2] [k/2]-l 

k=0 1=0 s=0 


n—4s 


2 n 


n + 2kjk —1\ I 


k f k — l\ f k — 2l\ f2s 
2 s 


(47) 


This is similar to fl38|l and fl43|l but with a different structure. Clearly, the fact that g = 4 
is not prime is related to the different appearance of the result. 

The case g = 5 presents some qualitatively new features. It is the hrst case for which 
the characteristic equation involves non-rational coefficients. Further, since this equation 
is quintic, it has (in principle) no analytical solutions. We will present no explicit formula 
for Tr and leave its full treatment for a future publication. 

Finally, we deal with the simpler g = 6 case. The solutions to the characteristic 
equation are now 

Ap = ±2^1 V^cos +^2p7r ^ p = (48) 

with 

36 -|- Cg 

cos ffl = - ^ . 

32^/2 

A calculation along the lines of the previous ones gives 

k/2] [^] [1/2] [l/2]-s 

Tr = EEEE ^ (^n—k—Al+bs—At ‘^21—As—At 

k=0 1=1 s=0 t=0 

{ n \{2k + 3l\ _l_ n -s\{l- 2sA /2A ^ 

\2k -)-3Z/\ k ) I — s\ s / \ 2t / V^/ 

which is even more elaborate than the previous formulas. 

This kind of computation can in principle be extended to higher values of g. Overall, 
we see no discernible pattern in the formulas obtained so far for Tr However, the 

fact that multiple not manifestly equivalent expressions exist for each g leaves the hope 
that a form more amenable to generalization may still exist. 


(50) 



5 Conclusion 

To conclude, we have further explored the previously-considered generating function 
^mi,m 2 /i,/ 2 (Q) of fho probability distribution of the algebraic area of random walks on 
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a square lattice. We have: (i) elucidated the relation of (“staircase walks”) 

to the cyclic sieving phenomenon (by demonstrating that its values correspond to the 
numbers of subsets hxed by powers of the cyclic permutation); (ii) shown that the values 
of (“biased walks”) are integers under certain conditions, and explained 

why; (iii) obtained an explicit expression for (Stembridge case), and con¬ 

nected it with the expression for the n-th moment of the Hofstadter Hamiltonian at the 
simplest nontrivial flux, q = 2] (iv) found explicit expressions for the same u-th moment 
for q = 3,4,6. Obtaining either a general closed expression for 

least more particular cases thereof, as well as a general form of the n-th moment evaluated 
at any root of unity, are open tasks. Note that the values of Zjni,m 2 ,ii,i 2 i^'^™^'^) 
numerically to be not integer anymore when q > 5 (to the exception of g = 6), making 
the task even more challenging. 
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